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Abstract 

The dynamics of piiantom bead-spring chains with the topology of a symmetric star with / arms 
and tadpoles (/ = 3, a special case) is studied, in the overdamped limit. In the simplified case 
where the hydrodynamic radius of the central monomer is / times as heavy as the other beads, 
we determine their dynamical eigenmodes exactly, along the lines of the Rouse modes for linear 
bead-spring chains. These eigenmodes allow full analytical calculations of virtually any dynamical 
quantity. As examples we determine the radius of gyration, the mean square displacement of a 
tagged monomer, and, for star polymers, the autocorrelation function of the vector that spans from 
the center of the star to a bead on one of the arms. 

PACS numbers: 05.40.-a, 02.50.Ey, 36.20.-r, 82.35.Lr 
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I. INTRODUCTION 



Bead-spring models play a central role in the theory and modelling of polymer dynamics. 
Most of the applications of bead-spring models (and polymer dynamics in general) is found 
for linear polymers, for which the polymer consists of a linear sequence of beads connected 
by harmonic springs. For a linear bead-spring polymer chain, with position R„ of the n-th 
bead, n = . . . N, the potential energy is thus 

N 

2 



n=l 

in which k is the spring constant. Besides the heavily studied linear polymer chains there 
are many other types of polymers whose dynamics deserve closer inspection. In this paper 
we concern ourselves with the dynamics of bead-spring chains that have the topology of a 
symmetric star (with / arms) and tadpoles. 

Further, in the context of bead-spring models, the Rouse model for linear phantom chains 
{l| deserves a special mention. In the Rouse model the dynamics of the beads is formulated 
in the overdamped limit, with a solvent with viscous drag coefficient and with thermal 
forces gn on the n-th bead, such that the dynamical equations of motion become 

dR„ 1 dU 



dt C dRn ' 
Here, the thermal forces are delta-correlated, namely 



(2) 



(g.(t)-g,(t')) = ^M(^-0, (3) 

with Boltzmann constant ks and temperature T. 

The reason why, to date, the Rouse model deserves a special mention, lies in the fact 
that it allows full analytical calculations of virtually any dynamical quantity. Rather than 
the equations of motion for the individual beads, one considers the so-called Rouse modes, 
whose amplitudes at time t are given by 

^ ^7i{n + l/2)p' 



1 r 



n=0 L -I 

with p = . . . N. The dynamics of these modes are uncorrelated, and one can derive the 

m 



following relation for p,q ^ exactly 



N 

Xpg{t) = (Xp(t) ■ X,(0)) = Ai- exp 



from the equation of motion (jl]): 



2 

-A2^t 5pg, (5) 



where Ai = 3kBT/(27T'^k) and A2 = 7r'^k/(. Here, and all throughout this paper, the angular 
brackets represent an average over the equilibrium ensemble of polymer chains. Equation 
([5D is further supplemented by Xop{t) = for p 7^ 0, and Xoo{t) = ([Xo(t) - Xo(0)]2) = 
6kBTt/{(N), where Xo(t) is the location of the center-of-mass of the polymer at time t. 



Using these mode amplitude correlation functions, the quantities of interest 



br a phantom 



chain can be analytically tracked by reconstructing them from the modes 

Using the definition of Rouse modes it has been shown in recent works that the dynamics 
of a tagged bead in a linear bead-spring model is described by the Generalized Langevin 
Equation (GLE) , and that the dynamics of polymers with steric repulsion (also known 
as self-avoiding polymer chains) as well as of those in melts as described by the repton 
model [3| can be well-approximated. Here, we continue this line of research, but now we are 
interested in exact solutions of the dynamical properties of polymers with the topology of 
stars and tadpoles. 

The structure of this paper is as follows. In Sec. |TT]we present the dynamical eigenmodes 
of symmetric phantom star polymers, and use the mode amplitudes to provide analytical 
expressions to the radius of gyration, mean square displacement of a tagged monomer, and 
the autocorrelation function of the vector that spans from the center of the star to a bead 
on one of the arms. In Sec. IIIII we repeat the exercise for chains with tadpole topology. We 
end the paper with a short discussion in Sec. IIVI 



II. DYNAMICAL EIGENMODES OF SYMMETRIC STAR POLYMERS 

A major difficulty for dynamics of polymers with a more complex topology like the 
symmetric star polymer is that in most cases an elegant analytical expression for the set 
of dynamical eigenmodes cannot be found jsj. Here we show that for a symmetric star 
polymer with a special central bead the dynamical eigenmodes can indeed be written down 
exactly, and subsequently the dynamical behavior of many interesting physical quantities 
can be determined precisely. For simplicity, in commensuration with the Rouse model, 
henceforth we term these dynamical eigenmodes as Rouse modes. Specifically, we consider 
a star polymer with / identical arms, each consisting of N identical beads, connected to a 
central bead whose hydrodynamic radius is / times as big as the other ones, i.e., its viscous 
drag is / times as large. A graphical representation of such a star polymer with / = 5 and 
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= 4 can be found in Fig. [T] It is worthwhile to note in this context that a key method 
to synthesize a star polymer chain is to attach the arms, which are linear chains, to a 
multivalent central core with sticky ends (see, e.g., Ref. Q and the references cited therein). 
Although from the synthesis process it is realistic that the core has a significantly higher 
hydrodynamic radius, the choice in our simplified model to make hydrodynamic radius of 
the central bead exactly / times as the other beads is motivated by our strive to determine 
the Rouse modes exactly. For modestly-sized chains, where the hydrodynamic radius of 
the core is not / times as big as the arm beads, the dynamical matrix (the homogenous 
part of the dynamical differential equation) has been diagonalized numerically, yielding the 
numerical identification of the Rouse modes jsj. 



FIG. 1: Two graphical representations of a symmetrical star polymer to visualize the Rouse modes. 
Depicted is a symmetric starpolymer with five arms consisting of four beads each and a central 
monomer which has a friction coefficient five times that of beads in the arms. The transparent 
configurations are polymers in the origin stretched in the xy-plane for visual convenience with no 
Rouse mode excited in the z-direction. The opague configurations are like the transparent ones but 
with a pure Xi mode (a) and y[*'"'^ modes (b) in the z-direction. 

We label the position of the central bead as Rq and let Ra,„ be the position of the n-th 
bead, n = 1 . . . N, in the a-th arm, a = 1 . . . f . We consider two types of Rouse modes given 




(a) 
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by 



Xp(t) 



N- 



cos 



7rp/2 
N + 1 



1 ^ 

Ecos 

n=l 



^ f,N 
a,n=l 

7r(A^-n + l/2)(g- 1/2) 
AT + 1/2 



7r(n + l/2)p' 



R-a,n(^) 



+ 1 



(6a) 
(6b) 



with p = . . . N and q = 1 . . . N . A visualization of the two kinds of modes can also be 
found in Fig. [1] The first set of modes Xp(t) are like Rouse modes for a stringpolymer 
through all the arms. The second set of modes Yp'''\t) can also be thought of as Rouse 
modes as in Eq. (|4]) with p odd valued and through a linear chain of length 2A^ + 1 made up 
by arms i and j and the central bead. There are Nf{f — l)/2 modes of the type Yp'''\t) 
with i < j, but the total set of these modes contain only (/ — 1)A^ independent degrees of 



freedom, for instance because we could have constructed, for 

every p, an orthogonal set of / — 1 modes out of the full set of modes Yp'''\t), but we choose 
not to do that for the sake of mathematical elegance. Combined with + 1 modes Xp(t), 
the total set contains fN + 1 three-dimensional modes needed to describe the system that 
has fN + 1 beads so that the number of degrees of freedom is the same. 

Similarly to the single chain, the dynamics of these modes for a star polymer with long 
arms are captured by 



Xp,{t) ^ {Xp{t) ■ X,(0)) 



3hT N 



27r2A; /p2 
Y(M)(i).YfO(o)) 



exp 



5, 



3hT 



N 



exp 



7r^k{p-l/2f 



C 



(7a) 



(7b) 



SiT^k (p- 1/2)2 

where = Sik - Sjk - Su + 6ji. This is supplemented by Xoo(t) = ([Xo(t) - Xo(0)]2) = 

6kBTt/{(f{N+ 1)) and all other correlations between modes strictly zero. The derivations 
of the Rouse mode amplitudes can be found in Appendix A. 



A. Radius of gyration 

The squared radius of gyration is defined as the weighted sum over all differences between 
the position of a monomer and the center of mass squared. Below we work out the radius of 
gyration for the case of the star polymer where the central monomer is / times heavier than 
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the other beads, although the radius of gyration can be calculated following the same line 
for other well. In the former case the location of the center-of-mass Hcmit) = Xo(t), 

and the squared radius of gyration is defined as 



Ri 



1 



f{N + l) 



f,N 



i,n=l 



which can be calculated by plugging in Eq. (lAip . This reduces to 



N 



NJ,fJ 



p=l •' ^ ' p,i,j,k=l 



Plugging in Eq. ([7]) explicitly and taking the long polymer limit yields 

fJJ 

"9 ~ ' 



kir"^ f ^ — ' 
•' p=i ^ 



(8) 



(9) 



The sums can be evaluated by using Eq. (lB4p . 

And so the radius of gyration squared for a long symmetric star polymer becomes 

.^ kBTN2>f-2 
k 2f ■ 

This result is consistent with that of a linear chain (/ = 1). 



(10) 



(11) 



B. Mean square displacement of a monomer 

We now consider the mean square displacement of the central bead and a bead in an arm 
for a star polymer with long arms. We start by defining the displacement vector for the 
central bead and writing it in terms of modes using Eq. (lAlaP : 



ARo(t) =Ro(t) - Ro(0) 

=Xo(t) -Xo(0) + 2^cos 



N 



p=l 

The mean square displacement is then given by 

AT 

(AR2(t)) =Xoo(t) + 8j]cos2 

p=i 



np/2 
N +1 



{X,(t)-X,(0)} 



(12) 



7rp/2 



{X,,(0)-X,,(t)}, 



(13) 



where Eq. (17al) can be plugged in, and the orthogonality of the modes was already used 
for simplification. Before doing so, let us look at very short time scales t < (/{4k). The 
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exponent in Xpq{t) can then be expanded and the sum exactly evaluated using Eq. (lB2ap . 
resulting in (ARQ(t)) = QkBTt/{(f) dominating over the mean square displacement of the 
whole polymer which at short time scales is negligible. At very long time scales Xpg(t) goes 
to zero and the term with the summation in Eq. (1131) can again be exactly evaluated using 
Eq. (lB3ap . So for t < (N^/{3k) the summation has a larger contribution than the mean 
square displacement of the whole polymer. For intermediate times the summation dominates 
and can be rewritten as an integral for very long polymers: 



exp 



c 



-X 



UkeT ri 



(14) 



The mean square displacement as approximated for intermediate times is greater than that of 
the approximation for the short time scales for t > A^/lnk), so that is when the intermediate 
time regime begins. For the central bead in a star polymer the mean square displacement 
then becomes 



(ARg(t)) 



( QksT 
UkeT 



fort < 



4k 



t 



f V ^kC 
QkeT 



for — < t < — — 



(15) 



fort > 



~3k 



Figure [2] shows the exact evaluation of Eq. (fT3|) for some star polymer together with ap- 
proximations made for the short, intermediate and long time scales found in Eq. f lT5|) . The 
central bead first behaves like a single bead with friction coefficient Cf . After that the 
movement is restricted by local connections to surrounding beads in the polymer. For very 
long times the position of the bead within the polymer is negligible and the mean square 
displacement behaves as that of a single bead with friction coefficient Cf{N + 1). 

A very similar approach can be used for the mean square displacement of a bead in an 
arm of the star polymer, defined as 



AR„(t) =Ri,n{t) - Ri,„(0) 



N 



=Xo(t) -Xo(0) + 2^cos 
p=i 

f,N 

COS 



4 



i>p=i 



7r(n + l/2)p' 
N + 1 

7r(iV-n + l/2)(j9 - 1/2) 



iV+ 1/2 



{X,(t)-X,(0)} 
{Y(^'^)(t)-Y(^'^)(0)}. 



(16) 



Using orthogonality of mode Yp with X^ and with Yg^'^'' with q p, and evaluating the 




t 

FIG. 2: The scaled mean square displacement {AJi? (t)) / \/t as a fmiction of time. The mean 
square displacement of several beads in a symmetric star polymer as given by Eqs. (jl3ll7p were 
exactly evaluated with A = 12 and N = 10^ and other parameters put to 1. The solid blue line 
corresponds to the mean square displacement of the central bead whereas the red, green, magenta, 
and cyan solid lines correspond to n = 10"'^ A^, N/2, 9A^/10, and N in Eq. (jl7p respectively. A bead 
positioned somewhere along the arm will at first behave as if it were in the middle of the arm. After 
time grows either the end or the center will become an influence at first and the bead will start 
mimicking a bead at one of those places. The short time scale t < C,/{Ak) and the very long time 
scale t S> N'^Ci'^f ~ 2)/(3A;) for which (AR^(t)) ~ t corresponding to the dotted and dot-dashed 
lines are separated by a time during which (AR^(t)) ~ \/i corresponding to the dashed line as in 
agreement with Eqs. (jl5|19ti20|) . 
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double sum over the arms, the mean square displacement becomes 



N 



(AR2(t))=Xoo(t) + 8 5^cos^ 



16(/ 



/ 



p=i 



7r(n + l/2)p' 



cos 



p=l 



N + 1 

7i{N -n + l/2){p-l/2) 



N+1/2 



x|ya,2)(i,2)(o)_^am2)(,)|. 



(17) 



For small values of t the exponents can again be expanded and the sum exactly evaluated 
using Eqs. flB2atlB2b|) so that in the long polymer limit (AR^^(t)) = QksTt/C,. The sum- 
mation reaches its maximum, using Eqs. (lB3atlB3b|) . on long time scales so that for a bead 
at the end of an arm the mean square displacement term of the whole polymer will start 
to dominate for t > N'^({3f — 2) /{3k). The closer a bead is to the central bead the faster 
its mean square displacement will behave like that of the whole polymer. We can again 
approximate the summation by an integral for intermediate times but depending on the 
position of the bead in an arm it will behave differently. 



(AR^(t)) 



UkeT 
' -K^kf 



°° dx 

~~2 



1 — exp 



X < cos 



^' 

7r(n + l/2)p' 



knH 



C 



-X 



+ (/-!) cos 



7r(A^-n + l/2)(p- 1/2) 
iV+ 1/2 



(18) 



The exponent in the integral will suppress the contribution for higher p-values. For a bead 
at the end of an arm, n equals A^, the cosines can be taken to be 1 and the integral as the 
same as the integral for the central bead but with an extra factor /. For a bead in the 
middle of an arm, n equals A^/2, in the long polymer limit the first cosine will be 1 for even 
values of p and for odd values, and the second cosine will be 1/2 for all values of p. This 
will give the same integral as for the bead at the end of an arm but smaller by a factor 1/2. 
We thus find for the mean square displacement of a bead at the end of an arm 



c 



(AR^(t)) = <^ l2kBT 



fort < 



Ak 



GkeT 



TikC 



fori^<t<^ 
irk 3k 

fort > ^(3/ -2) 



(19) 
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For the bead located exactly in the middle of an arm the mean square displacement becomes 



(AR^/2(t)) 



( eksT 
GkeT 



c 

t, for t < -V 
Ak 

t , AC N\ 

— — , f or — < t < — — 

TikC, nk 3k 

t, fort > ^(3/ -2) 



(20) 



Figure [2] shows the exact evaluation of Eq. f lT7|) and demonstrates the validity of Eqs. fll9t - 
20|) . The behavior of a bead somewhere along the arm can also be extracted from these 
results. At first it will behave like the bead in the middle of an arm since locally they are 
the same. As time progresses it will either start feeling the end or the center of the polymer 
at first and mimic the behavior of the bead at the end or center respectively after which the 
mean square displacement will behave like that of the whole polymer. 



C. Correlation function of a vector connecting a bead to the central bead 



Consider the spatial vector connecting a bead in some arm to the central bead 



The correlation function C„(t) = (rj^„(t) ■ rj„(0)) is then given by 



(21) 



N 



C„(t)=16^si 



sm 



p=i 



Ti{n + l)p/2 
N +1 



sm 



Tinp/2 
N +1 



Xpp{t) 



16 



cos 



7r{N -n + l/2){p- 1/2)' 
1/2 



(22) 



j,k,p=l 

where the cosines were reduced to sines for notational convenience. Having filled in the 
mode dynamics functions explicitly yields 



2AkBT 



N 



Cf{N 



sm 



p=i 



7r(n + l)p/2 



N + 1 



sm 



Tmp/2 
N +1 



■ exp [-axpt] 



+ 



12(/-l)fcgT 
C/(2iV+l) 



N 



cos 



p=l 



7i{N -n + l/2){p- 1/2) 
1/2 



exp [-aYpt] , (23) 



where and ctYp are defined in Eqs. (lA5IIA6p . At very large and very small times the 
exact correlation functions for the modes in Eq. flASp can be used. At very small times the 
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10-4 

10-2 10-1 10^ 10^ 10^ 10^ 10^ 10^ 10^ 10^ 

t 

FIG. 3: The scaled correlation function Cn{t)/n as defined in Eqs. (|2m22p in a double-logarithmic 
plot as a function of time. The correlation functions for several monomers in a symmetric star 
polymer as given by Eq. (|23p were exactly evaluated with / = 12 arms of length N = 10^ and 
other parameters put to 1. The red disks, green triangles, and blue squares correspond to n = 1, 10, 
and N for Cn{t) respectively. Approximations for the correlation function in the short time scale 
l)^/(4vr2A;), intermediate time scale Cn^{f - lf/{kTrf^) < t < N^C/{Tr^k), and the 
very long time scale t > iV^^/(7r^/c) as in Eq. (|25p correspond solid lines. The approximation for 
intermediate times is only very accurate for small n, but for larger n this time domain becomes 
smaller or even non-existent as is the case for n equals as can be seen in the figure. 



exponential can be omitted and the resulting equations solved using Eqs. (lB3b|B3dp . Adding 



the two terms then results in C„(t) = 3nkBT/k. To determine the region for which this 
approximation is valid we notice that the exponent is roughly 1 for < ( N"^ / {tt'^ kt) . The 
term in the summation will go to for p = 2{N -|- l)/(?7, -|- 1) and so the largest contribution 
to the summation is for all the terms before this happens. Solving this for t where the large 
polymer limit is taken gives t < ([n + l)^/(47r^A;) for which the approximation is valid. For 
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very large time values t > {N'^(){7r'^k) only the lowest mode will give contribution. Since 
/ = 1,2 are two cases of linear chains we focus on / > 3 for which the second term in 
Eq. (123|) will dominate. Taking the long polymer limit the sine can be expanded and the 
cosine reduced for notational convenience. For very large times the correlation function can 
thus be approximated by = l2kBTN/{Tx^k) (1 - cos [7m/A^])exp [-Ti^kt/ {ACN^)]. For 

intermediate time values the second term can be approximated by a Gaussian integral in 
the long polymer limit by expanding the cosine for small n. 



6{f-l)n^kBT 
kf 



dx exp 



C 



-X 



3(/ 



(24) 



For t > (n'^if — lY / {nkp) the intermediate time approximation will be smaller and thus 
more accurate than the approximation for small times. Putting this together gives for / > 3 



Cn{t) 



SnksT 

3{f-l)n'^kBT fkt 

kJT^ U 
l2kBTN 



-1/2 



for t < — -— — 
ATT^k 

(n^if - 1)2 



TT 



'■k 



COS 



7rn 



N 



exp 



TT^k 



for ■ 



fort > 



kTTp 



(25) 



A graphical representation of the correlation function for a spatial vector between some 
monomer and the central monomer in a symmetric starpolymer with / = 12 and = 10^ 
is shown in Fig. |3l 



III. DYNAMICAL EIGENMODES OF TADPOLE POLYMERS 

A tadpole polymer can be seen as a star polymer where two of the three arms have the 
ends connected. Here we extend our calculations of the previous section and write down 
the exact solution for the eigenmodes for the specific case where the tadpole is built from 
a symmetric star polymer with arms of length A^ and a central bead with hydrodynamic 
radius three times as large as that of all the other beads. The Rouse modes are then a 
variation of the modes for a ring polymer combined with those for a star polymer. The first 
set of A^ + 1 modes are like the X modes for the star polymer where all the arms behave the 
same. The second set of A^ modes are very similar to the Y modes of the star polymer where 
the ring takes on the role of an arm and the tail of the tadpole the role of another arm. The 
third set Z of A^ modes are Rouse modes for a ring polymer but specifically such that it is 
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antisymmetric around the central bead. In Fig. H] these three different modes are depicted 
for p = 1. We label the tadpole polymer like a three-armed symmetric star polymer where 




FIG. 4: Three graphical representations of a tadpole polymer to visualize the Rouse modes. De- 
picted is a tadpole polymer, which can be seen as a symmetric star polymer with three arms 
consisting of four beads each where two arms are connected at the ends, with a central monomer 
which has a friction coefficient three times that of beads in the arms. The transparent configu- 
rations are polymers in the origin stretched in the xy-plane for visual convenience with no Rouse 
mode excited in the z-direction. The opague configurations are like the transparent ones but with 
a pure Xi mode (a), Yi mode (b), and Zi mode (c) in the z-direction. 



arm one and two have the ends connected. The Rouse modes are then given by 

3,N 

3 



X,(t) 

z,(t) 



1 


N + 




1 




AN- 




1 





cos 



np/2 



N+1 



2N + 



N 

n=l 

N 



COS 



(t) + I COS 
a,n=l 

7r(iV-n + l/2)(g-l/2) 



7i{n + l/2)p 



(26a) 



sm 



n=l 



1/2 

7r(iV-n + l/2)g' 
N+1/2 



N + 1 

Rl,n(t) + R2,n(t) - 2R3,„(t)],(26b) 
Rl,n(t) - R2,n(t)) , (26c) 
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with p = . . . N and q = 1 . . . N. The vahdity of these Rouse modes as the dynamical 
eigenmodes can be proved hke for the star polymer in Appendix A. Let us define 

N 



X,(t) = Xo(t) + 2 J]cos 
p=i 



7r{n + l/2)p 



N 



%=1 



N 



p=l 



N + 1 

7r(iV-n + l/2)(p - 1/2) 
iV+ 1/2 

7r{N -n + l/2)p 



Yp(t), 



iV + 1/2 



Zp(t), 



(27a) 
(27b) 
(27c) 



so that the bead locations are given by 

Ro(t) 

R-l,n.(^) 
R-2,ri(i) 
R'3,n(^) 



Xo(t) 

X„(t)+Y„(t) + Z,(t), 
X„(t)+Y„(t) -Z„(t), 
X„(t)-2Y„(t). 



(28a) 
(28b) 
(28c) 
(28d) 



The dynamics for these modes are very similar as for the modes of the star polymer. The 
non-vanishing correlation functions for the modes of the tadpole are then given by 



([Xo(t)-Xo(0)]2) 
(X,(t) ■ X,(0)) 
(Yp(t)-Y,(0)) 



t 



2kBT 
C{N + 1) 
ksT 1 

C(Ar + l)2ax, 
dkeT 1 

4C(2A^ + l)2aY, 



exp [-axp^] Spg 
exp [-aYpt] Spq 



(29a) 
(29b) 
(29c) 
(29d) 



for p,q = 1, . . . N and 



axp =4- sin^ 



2iV + 2 



, "Yp = 4^ sin^ 



7r(p-l/2) 
2A^ + 1 



. ^ . 2 

«Zp = 4- sm 



Trp 



2N + 1 



(30) 



A. Radius of gyration 

When the ends of two arms of a three armed symmetric star polymer are connected the 
mobility reduces and the radius of gyration will become smaller. Following the same steps 
as for the star polymer the radius of gyration is defined in the same way for the case when 
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the central bead is three times as heavy as any other bead, and by plugging in the inverses 
it can be rewritten in terms of summations over the correlation functions of the modes at 
time zero. The radius of gyration squared for a long tadpole becomes: 

- ^^^^ 
which is 5/7th of the radius of gyration of the same polymer but with the ends not connected. 



B. Mean square displacement of a monomer 

Since the tadpole can be seen as a three armed star polymer with ends of two of the arms 
connected, the mean square displacement of the monomers in the third arm and that of the 
central monomer will not feel the effect of the connection of the two remaining arms. And 
so for the central monomer and the monomers in the tail of the tadpole the expressions will 
be exactly the same as in Eqs. (I13|17p respectively. A monomer in the closed circle of the 
tadpole should behave in the same way as some monomer roughly in the middle of the tail: 



ARo(t) = Ro(t) - Ro(0) , AR,,„(t) = R,,„(t) - R.,„(0). 



(32) 



The mean square displacement for monomers in a tadpole can be described by 
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FIG. 5: The scaled mean square displacement (AR^(t)) /\/t as a function of time. The mean square 
displacement of several beads in a tadpole made by a symmetric star polymer with ends of arm 1 and 
2 connected as defined by Eq. ([32|) were exactly evaluated with N = 10"* and other parameters put 
to 1. The solid blue line corresponds to the mean square displacement of the central bead whereas 
the red, green, magenta, and cyan solid lines correspond to i = 1 and n = lO^^N, N/2,9N/10, 
and N in Eq. (f32|) respectively. Unlike for the unconnected arms in the star polymer here there 
is no real difference between a monomer at the middle or end of the connected arm. A bead 
positioned somewhere along the arm will at first behave as if it were in the middle of the arm. If 
it is close enough to the heavy central bead the movement will be restricted and the mean square 
displacement will mimic that of the central bead. The short time scale t < (^/(Ak) and the very 
long time scale t ^ 5N'^C/{6k) for which (AR^(t)) ~ t corresponding to the dotted and dot-dashed 
lines are separated by a time during which (AR^(t)) ~ \/t corresponding to the dashed line as in 
agreement with Eqs. (j33p . 
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Note that the domains of vahdity for the different behavior have changed in comparison to 
the star polymer. The exact sums for the monomers in a tadpole have a slightly different 
maximum as can be calculated using Eq. (IB3p . 

IV. DISCUSSION 

For bead-spring models of polymers with the topology of a symmetric star with / arms, 
where the hydrodynamic radius of the central bead is / times as heavy as any other bead, we 
derived the exact expressions for the dynamical eigenmodes. We demonstrated the usefulness 
of this exercise by exact calculations of various quantities that yield prefactors as well as 
give insights into the scaling behavior at different time regimes for these quantities. The 
radius of gyration was shown to be given by = ksTNi^if — 2)/{2kf); the mean square 
displacement of the central bead and of various other individual beads was shown to scale 
as ~ t at very short as well as at very long times, with an intermediate regime in which it 
scales as ~ \/t; and the correlation function of an orientational vector was shown to stay 
invariant at very short times, decay exponentially at very long times, and decay as ~ 1/ \/i 
in the intermediate regime. Similar results were also derived for tadpoles, which are / = 3 
star polymers in which two arms are connected. 

Our exact results are limited to the very specific case of /-arm star polymers and tadpoles 
in which the hydrodynamic radius of the central bead is exactly / times bigger than the other 
beads. The Rouse modes can still be constructed for other cases, since the Rouse equations 
are linear in the bead positions; indeed, this has already been achieved numerically^]. We 
do not rule out analytical forms of the Rouse modes for these cases, although we can expect 
that these analytical forms may not look aesthetically pleasant. Even if the hydrodynamic 
radius of the central bead is not exactly / times bigger than the other beads, we do not forsee 
differences in the behavior of the mean-square displacements in the scaling limit, namely 
the three regimes mentioned above, albeit with different prefactors. The same should hold 
for modest modifications of the interatomic bonds: no qualitatively different behavior is 
expected, but prefactors will be affected. As for the dynamical eigenmodes of star and 
tadpole polymers, qualitative differences will however arise due to significant hydrodynamic 
interactions between the beads, excluded- volume effects, or if the solution is no longer dilute. 

Apart from these direct results, the usefulness of our exercise lies perhaps even more 
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in providing a proper language in which to characterize the dynamics of star and tadpole 
polymers; very similar to the ubiquitous use of Rouse mode amplitudes for linear polymers. 

In future work, we will use the mathematical expressions for the Rouse modes derived 
here, to characterize the dynamics of star polymers in simulations of self-avoiding stars, 
in dilute circumstances, as well as in various dense environments (such as a melt of linear 
polymers, or a homodisperse melt of star polymers). 



Appendix A: Determination of the Rouse modes of a star polymer 



Several steps are needed to derive the dynamics of the modes ([7]) from the equations of 
motion for a symmetric star polymer. It is useful to note that the modes are complete, 
which allows us to express the positions of the beads from the mode amplitudes, given by 
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The correctness of these equations can be checked using the orthogonality relations (IBip . 

The dynamical equations of motion for the star polymer combined with the potential 
energy then result in the following equation of motion for the beads: 
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Note that i = 1 . . . / and Eq. (lA2cP is valid for all n = 2 ... where Eq. (lAlbp is needed 
to show that Rj^Ar+i = Rj^at. For convenience we define 
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(A3) 



for n = Q . . . N which coincides with the inverse of Rq as given in Eq. (lATa]) . Taking the 
time derivative on both sides of Eq. (l6a|l . plugging in the equations of motion for the beads 
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from Eq. (lA2p . and using the definition above for which R_i 
in 
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where Gp is the transform of the thermal forces which will be calculated later. By using (lA3p . 
the trigonometric identities, namely the angle sum and difference identities and the power- 
reduction formula, and finally the orthogonality relation the set of differential equations 
becomes 
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The set of differential equations for the Yp modes can be written down in a similar manner. 
We define Rn'"''' = Rj_„ — Rj^„ and follow similar steps as for the Xp modes so that 
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for some transform of the thermal forces Gp'"'^ 



With the above transformations, the set of differential equations for the beads have been 
transformed to a set of disconnected linear differential equations. 

To find the relation for the modes as in Eq. ([7]) we must first determine the transform of 
the thermal forces Gp and Gp'"'^ The transforms are equal to that of the modes in Eq. ([6]) 
but by replacing the position of the beads by their thermal force. Use that the thermal forces 
are uncorrelated in time and between beads as in Eq. ([3]) and recall that the central bead 
has a friction coefficient / times as large as that of the other beads. The only nonvanishing 
functions with p = 1 . . . N are 
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With these, the set of differential equations can be solved exactly resulting in the following 
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relations between modes: 
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where all other correlations between modes are strictly zero. By taking the long-polymer 
limit the sines can be expanded up to second order and the results are in Eq. ([7]). 



Appendix B: Useful mathematical relations 



The following relations are useful for calculating the Rouse modes: 
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for p,q = 1 . . .N and where the left hand side of Eq. (IB lap equals 26po for p = . . . N. 
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